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Orthogonal polynomial random matrix models of N x N hermitian matrices lead to Fredholm 
determinants of integral operators with kernel of the form {(fi{x)^p{y) — tp{x)ifi{y)) /x — y. This paper 
is concerned with the Fredholm determinants of integral operators having kernel of this form and 
where the underlying set is the union of intervals J = IJJL]^ (12^-1, a2i). The emphasis is on the 
determinants thought of as functions of the end-points at- 

We show that these Fredholm determinants with kernels of the general form described above are 
expressible in terms of solutions of systems of PDE's as long as (p and satisfy a certain type of 
differentiation formula. The {ip, tp) pairs for the sine, Airy, and Bessel kernels satisfy such relations, 
as do the pairs which arise in the finite N Hermite, Laguerre and Jacobi ensembles and in matrix 
i-P ' models of 2D quantum gravity. Therefore we shall be able to write down the systems of PDE's for 

these ensembles as special cases of the general system. 

An analysis of these equations will lead to explicit representations in terms of Painleve transcen- 
CNJ ' dents for the distribution functions of the largest and smallest eigenvalues in the finite N Hermite 

and Laguerre ensembles, and for the distribution functions of the largest and smallest singular 
values of rectangular matrices (of arbitrary dimensions) whose entries are independent identically 
^ ' distributed complex Gaussian variables. 

^ . There is also an exponential variant of the kernel in which the denominator is replaced by e^^ — e^^ , 

' where h is an arbitrary complex number. We shall find an analogous system of difi'erential equations 

\ in this setting. If & = i then we can interpret our operator as acting on (a subset of) the unit circle 

■ in the complex plane. As an application of this we shall write down a system of PDE's for Dyson's 
' circular ensemble of N x N unitary matrices, and then an ODE if J is an arc of the circle. 

; I. INTRODUCTION 

D , It is a fundamental result of Gaudin and Mehta that orthogonal polynomial random matrix models of N x N 

■ hermitian matrices lead to integral operators whose Fredholm determinants describe the statistics of the spacing of 



X 



eigenvalues |28 3^. Precisely, if a weight fimction is given, denote by {pk{x)} the sequence of polynomials 



orthonormal with respect to w{x) and set 



Then E(n; J), the probability that a matrix from the ensemble associated with w{x) has precisely n eigenvalues in 
the set J (n = 0, 1, • • • ) , is given by the formula 

E{n;J)^^— ^—det {I ~XK^)\ (1.1) 
where Kjq is the integral operator on J with kernel 

N-l 

KN{x,y) := ^ ipk{x)ipkiy). 

fe=0 

It follows from the ChristofFel-Darboux formula (c.f. (^^) below) that XK^ix,!/) is a particular case of a kernel of 
the general form 
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This paper is concerned with the Fredholm determinants of integral operators having kernel of this form and where 
the underlying set is the union of intervals 

m 

J ■■= |J(a2j-i,a2j). 

The emphasis is on the determinants thought of as functions of the end-points a^. 
If we denote the operator itself by K then it is well known that 

■^logdetiI-K) = i-l)''-^Riak,ak) (fc = l,---,2m) (1.3) 
oak 

where R{x, y), the resolvent kernel, is the kernel of K{I — K)^^ . This requires only that A = 1 not be an eigenvalue 
of K and that K{x, y) be smooth. Jimbo, Miwa, Mori and Sato showed for the "sine kernel" 

sin(a; — y) 
x-y 

that if we define 

rk,± ■■= (/-i^)-'e±"(afc) 

then the i?(afc, Ofc) are expressible in terms of the rk^± and that these in turn, as functions of the ai, • • • , a2m, satisfy 
a completely integrable system of partial differential equations. They deduced from this that in the special case 
when J is an interval of length s the logarithmic derivative with respect to s of the Fredholm determinant satisfied 
a Painleve differential equation. (More precisely, s times this logarithmic derivative satisfied the so-called a form of 
Py of Jimbo-Miwa-Okamoto p^ , p3[ .) We refer the reader to for a derivation of these results in the spirit of the 
present paper. The discovery that Painleve transcendents can be used to represent correlation functions in statistical 
mechanical models first appeared in the 2D Ising model |^,|2^,^ . 

The sine kernel arises by taking a scaling limit as — > oo in the bulk of the spectrum in a variety of random 
matrix models oi N x N hermitian matrices. But if we take the Gaussian unitary ensemble (also called the Hermite 
ensemble; see below), and others as well, and scale at the edge of the spectrum then we are led similarly to the "Airy 
kernel" 

Ai(x)Ai^(j/) - Ai'{x)Ai{y) 
x~y 



where Ai(a;) is the Airy function ||6|,^ 30 3^]. For this kernel the authors found [g8[ a completely analogous, although 
somewhat more complicated, system of PDE's, and showed that for J a semi-infinite interval (s, oo) there was also 
a Painleve equation associated with the determinant — this time P//. Similarly, if we scale the Laguerre ensemble at 
the left edge of the spectrum or the Jacobi ensemble a t ei ther edge (see below for these ensembles also), then we 



obtain yet another kernel, the "Bessel kernel," where in (1.2) (p{x) = Ja{y/x), ip{x) = xip'{x) with Ja the usual Bessel 



function. Again we found [g9[ a system of PDE's for general J and, for J = (0, s), a Painleve equation associated 
with the Fredholm determinant — this time Py (actually a special case of Py which is reducible to Pm). 

In looking for (and finding) analogous systems of PDE's for finite N matrix ensembles we realized that all we needed 
were differentiation formulas of a certain form for ip and ip, namely 

m{x)(p'{x) = A{x)ip{x) + B{x)'ip{x), 

m{x)4)'{x) = -C{x)ipix) - A{x)^p{x) (1.4) 

where m. A, B and C are polynomials. 

The ((yC, ij}) pairs for the sine. Airy, and Bessel kernels satisfy such relations {m{x) = 1 for sine and Airy, m{x) = x 
for Bessel) as do the pairs which arise in the finite N Hermite, Laguerre and Jacobi ensembles {m{x) = 1 for Hermite, 
m{x) = X for Laguerre and m{x) = 1 — x^ for Jacobi) and therefore we shall be able to write down the systems 
of PDE's for these ensembles at once as special cases of the general system. An analysis of these equations will 
lead in the cases of the finite N Hermite and Laguerre ensembles to explicit representations in terms of Painleve 
transcendents for the distribution functions for the largest and smallest eigenvalue. A consequence of the latter is 



2 



such a representation for the distribution functions of the largest and smallest singular values of rectangular matrices 
(of arbitrary dimensions) whose entries are independent identically distributed complex Gaussian variables; for these 
singular values are the eigenvalues of a matrix from an app ropriate Laguerre ensemble JTzt . 

There is also an exponential variant of the kernel (^]^) in which the denominator is replaced by e^^ — (or 
equivalently sinh|(a; — y)) where b is an arbitrary complex number. With an appropriate modification of ( |l.4| ) we 
shall find a completely analogous system of differntial equations. Observe that if 6 = i then we can interpret our 
operator as acting on (a subset of) the unit circle in the complex plane. As an application of this we shall write down 
a system of PDE's for Dyson's circular ensemble oi N x N unitary matrices, and then an ODE if J is an arc of the 
circle. In case b is purely real our results have application to the so-called g-Hermite ensemble Pj3l|. 

Here, now, is a more detailed description of the contents of the paper. 



A. The Differential Equations 

In Sec. we derive our general system of partial differential equations. To describe these equations we first define 
the functions 

Q := {I - K)-^ ip, P := {I - Ky^ t^, (1.5) 

which depend also, of course, on the parameters a^, and then 

qk--=Q{ak), Pk-^P{ak) (fc = 1, • • • , 2to) (1.6) 
Ui := (Q(x),a;V(a;)) , := ((5(a;), xXa;)) , := {P{x),x'ij{x)) (i = 0, 1, • • •) 

where the inner products are taken over J. T hese are the unknown functions in our system of PDE's. We shall see 
that for any operator with kernel of the form (L2) we have for the resolvent kernel the formulas 

i?(„,„,)^^Z^i^» (j^fc), (1.7a) 

R{ak,akj=Pk^ qk-^ — , (l-7b) 

oak oak 

for the Qj and pj the differentiation formulas 



^^{-iyR{a„ak)qk {j ^ k), (1.8) 
oak 

^^{-l)^R{a,,ak)pk ij^k), (1.9) 



and for the u.-, i;.-, Wj differentiation formulas of the form 



—-^ , — , —-^ = polynomial in pk , qk and the various Ui,Vi, wt 
Oak oak oak 



These equations are universal for any kernel of the form (1.2). What depends on (1.4) are the remaining differential 
formulas 

dq ■ 

m{aj)—-^ = polynomial in qj, pj and the Ui, Vi, Wi 



daj 



^(-l)'=i?(aj,afc)gfc, 



Op 

— polynomial in g,, n„ and the Ui, Vi, Wi 

da 



Y^{-lfR(a.j,ak)pk: 



and the representation 
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m(aj)R{aj , Gj) = polynomial in qj, pj and the Ui, Vi, Wi. 

The polynomials on the right sides are expressed in terms of the coefficients of the polynomials to, A, B, C in (|l.4|). 
We mention that in pst no "extra" quantities Ui, Vi, Wi appear, but this is quite special. In general the number of 
triples (ui ^Vi,Wi) which occur is at most 

max (deg A, deg B, deg C, deg to — 1) , 

although in practice fewer of these quantities actually appear. 



B. The Examples 



First in Sec. [I] we quickly derive, as special cases, the systems of equations for the sine. Airy and Bessel kernels. 
Then in Sec. [V we derive and investigate the equations for kernels "beyond Airy." To explain this we replace the 
variables x, y in (1^) by A and /z, think of (a completely new variable) x as a parameter, and observe that for each x 



ki{x + A)Ai'(x + /i) - Ai'(a; + A)Ai(a; + ^) 

A — yU 



(1.10) 



has the same properties as the Airy kernel. (In the differentiation formulas ( |l.4[ ) the variable is now A and a; is a 
parameter in the coefficients.) Observe also that Ai(a; + A) is, as a function of x, an eigenfunction of the Schrodinger 
operator with potential £,{x) = —x corresponding to eigenvalue A. 

In the hermitian matrix models of 2D quantum gravity [p|,[7| jl^Jl^ , ^ solutions to the so-called string equation 

[Q,V\ = 1 

determine the functions Lp and "0. In the simplest case of the KdV hierarchy, the operator Q is the Schrodinger 
operator (note our convention of sign of Dl) 



Dt 



and the differential operator V (in x) is 



V = 



|^g(2£-l)/2^ (^=1,2,...) 



where (•)+ is the differential operator part. The potential ^ then satisfies a differential equation determined by the 
string equation and ip{\, x) is the eigenfunction 



satisfying 



Setting 



the kernel 



is then 



{Dl+^{x))^{\x)^M\x) 



ip{\x) ^ Dxip{\x) 

ip{\, x)ip{n, x) - V'(A, x)(p{fi, x) 



(1.11) 



(1.12) 



(1.13) 



^(A,m) = 



A — 

Lp{X,y)(p{fi,y)dy . 



(1.14) 



These are the kernels which we say a re "b eyond Airy" since for £ = 1, V — D^, ^{x) = ~x, (1.14) reduces to the 
generalized Airy kernel ( 1.10 ). From ( 1.11 ) and ( 1.12 ) it follows that for general £ the hmctions (p{X,x) and '0(A,x) 
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satisfy difFerentiation formulas (in A) of the form ( |l.4| ). (Again in the differentiation formulas (1.4) the variable is 
now A and a: is a parameter in the coefficients.) In Sec. |V| we illustrate these methods for the case 1 = 1. 

In Sec. we study in some detail the finite A'' Hermite, Laguerre, Jacobi, and circular ensembles. In orthogonal 
polynomial ensembles one is given a weight function w{x) and then, for any symmetric function / on R^, we have 

E{f{Xi, ■ ■ ■ ,Xn)) = Cn J " J /(^l' ■ ■ • ,2:Ar) ^(Xi) Y[ - Xjf dxi ■ ■ -dxN (1-15) 

i<j 

where "£"' denotes the expected value, Ai, • • • , Ajv are the eigenvalues, and cat is a constant such that the right side 
equals one when / = 1. In the Hermite ensemble w{x) = and the integrations are over R, in the Laguerre 

ensemble w{x) = x°'e~^ and the integrations are over R+, and in the Jacobi ensemble w{x) = (1 — a;)"(l + x)^ and 
the integrations are over (—1, 1). In the circular ensemble w{x) ~ 1 and the integrations are over the unit circle. 

The size parameter N will appear only as a coefhcient parameter in the equations we obtain; and we find that the 
equations for both bulk and edge scaling limits emerge as limiting cases. Our equations also make the study of large 
N corrections to the scaling limits tractable. 

For the Hermite, Laguerre and Jacobi ensembles there are natural intervals depending upon a single parameter s — 
for Hermite J = (s, oo) or (— oo, s), for Laguerre J = (s, oo) or (0, s), and for Jacobi J = (s, 1) or J = (—1, s) — and 
in all these cases we shall find an associated ord inar y differential equation. For Hermite and Laguerre these will be 



of Painleve type. Observe that taking n = in (1.1) shows that the Fredholm determinant in each of these cases is 



precisely the distribution function for the largest eigenvalue, or 1 minus the distribution for the smallest eigenvalue. 



C. General Matrix Ensembles 

In this final section of the paper we show that there are differentiation formulas of the form ( |1.4| ) when Hermite, 
Laguerre, or Jacobi weights are multiplied by e"^*-^-* where V{x) is an arbitrary polynomial. (Of course it must 
be of such a form that the resulting integrals are convergent.) It is the finite N matrix models corresponding to 
certain V{x) which, in an appropriate double scaling limit at the edge of the spectrum, lead to the kernels beyond 
Airy. (Strictly speaking this is true only for the universality classes ^ = 1, 3, 5, . . . as it is well-known that the cases 

= 2, 4, 6, . . . require coefficients in V{x) that make e"^'-^^ unbounded.) 

II. THE GENERAL SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS 

In this section we derive the system of partial differential equations that determine the functional dependence of 



the Fredholm determinant det( J — K) upo n th e parameters where K has kernel (1.2). After some preliminary 



definitions and identities in Sec. [I A, in Sec. pB| we derive those equations which are independent of the differentiation 



formulas (1.4). In Sec. II C we assume Lp and satisfy the differentiation formulas for the case m(x) = 1. Then in 



Sec. IIP we indicate the modifications necessary for the general case of polynomial m. Finally, in Sec. HE we derive 



the exponential variant of the system of equations. 

A. Preliminaries 

Our derivation will use, several times, the commutator identity 

[L,(/-i^)-i] = {I-K)-^[L,K]{I-K)-\ (2.1) 
which holds for arbitrary operators K and L, and the differentiation formula 

-f (/ _ K)-' = [I- K)-' ^ (/ - K)-\ (2.2) 
da da 

which holds for an arbitrary operator depending smoothly on a parameter a. 

It will be convenient to think of our operator K as acting, not on J, but on a larger natural domain V and to have 
kernel 

K{x,y)Xj{y) (2.3) 
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where Xj is the characteristic function of J and K{x,y) is the kernel (1.2). For example, for the sine and Airy kernel 
D = R, for the Bessel kernel 2? = R'^, and for the Jacobi kernel V — (—1, 1). The set J will be a subset of V. We will 
continue to denote the resolvent kernel of K by R{x, y) and note that it is smooth in x but discontinuous at y — au- 
We will also need the distributional kernel 

p{x,y) = S{x -y) + R{x,y) 
of (/ — K) ^ . The quantities R{aj, ak) in Sec. I A are interpreted to mean 



lim R{aj,y), 
ye. J 

and similarly for pj and qj . 

The definitions of Ui, etc. must be modified. Before doing this it will be convenient to introduce 

Qj{x) (/ - Ky^ x^ip{x), Pj{x) := [I - K)'^ x^ij{x) 

which for j = reduce to ( pTsI ) (Qo — Q^ Pq — P)- Then we define 



(2.4) 



Uj := {Q.x^ipXj) 




(2.5a) 


Vj := {Q.x^^jXj) 




(2.5b) 


ij := {P,x^^Xj) 




(2.5c) 


Wj := {P.x^^Xj) 




(2.5d) 


where the inner product (•, •) is now over the domain V 


That these definitions for Uj, Vj and Wj 


are equivalent to 



those of Sec. [A and the above equalities are left as exercises. They follow from the fact that 



(/ - fx, ^ (/ - K)-^ fx, = { ^^'-/jV ^ ' 



(2.6) 



for any smooth /. Here K* is the transpose of the operator K. (Note that K takes smooth functions to smooth 
functions while its transpose takes distributions to distributions.) We also observe that vq — vq. A final bit of 
preliminary notation is 

L = L{x,y) 

which means the operator L has kernel L{x,y). 

B. The Universal Equations 

In this subsection M denotes multiplication by the independent variable. If we consider the commutator of M with 
K and use the representation (2.3), we have immediately 



and so by (2.1) 



[M,K\ = {ip{x)ij{y) - ij{x)ip{y))Xj{y), 
M, (/ - K)-^] = Q{x) {I - i;Xj{y) ~ P{x) (/ - K'Y^ ^Xj{y). 



(2.7) 



(The transpose here arises from the general fact that if i = U{x)V{y) then T1LT2 = TiU {x) T^V {y) .) It follows 
immediately that 

A/, (/ - K)~^ = [x- y)p{x, y) = {x- y)R{x, y), 
and hence referring to (2.6), we deduce 



Rix,y)^9Mm^^mM i^^^yejy 



x-y 
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In particular we have deduced (1.7a) (recall definitions ( |l.6D ) and the representation 

i?(x, x) = Q'{x) P{x) ~ P'{x) Q{x) {x e J). 



(2.8) 



We remark that the generality of this expression for R{x, y) was first, as far as the authors are aware, stressed by Its, 
et. al. though it appears, of course, in the context of the sine kernel in the earlier work of JMMS |p5[| . 
We have the easy fact that 



_9_ 



K = (-l)''K{x,ak)5{y-ak) 



and so by (gJ) 



— {I-K)-' = {-lfR[x,ak)p{y.au). 
oak 



(2.9) 



At this point we use the notations Q(a;, a), P{x^ a) for P{x ) an d Q{x \ r espectively, to remind ourselves that they are 
functions of a as well as x. We deduce immediately from ( ^.9|) and (1^) that 



d d 
-— g(a;, a) = {-lfR{x, ak)qk, i^P{x, a) = {~lfR{x, ak)pk- 
oak oak 



(2.10) 



Since qj — Q{aj, a) and pj — P{aj, a) this gives 



I* = (-l)'=i?(a„a,)g,, ^ = {-l)'^R{a,,ak)pk, {] ^ k). 
oak oak 



These are equations (1.8) and (1.9). We record for use below 



dq.j _ d ^ d 

daj dx daj ' 



dpi 



d d 



(2.11) 



To obtain (1.7b) observe that (2.8) gives 



R{ak,ak) = TT- [Qix, a)pk - P{x, a)qk] 
ox 



But the expression in brackets above vanishes identically when x = ak and so the above is equal to 



d 

-7^ — [Q{x, a)pk - P{x, a)qk\ 
oak 



If we use ( ^.10 ) in the computation of this partial derivative, ( |l.7b| ) results. 

We now show that the Vk can be expressed in terms of the other quantities Ui, Vi and Wi (we could do Vk just as 
well) and that the Pk and Qk can be expressed in terms of these quantities and P, Q. From 



X^ _ y'^ 

x-y 



J2 ^^y' 



i+j=k-i 

i,j>0 



we get 



i+j=k-i 

ij>0 



and so 



M'', (/ - KY 



{q.{x){I-K') 'y^^P{y)Xjiy)-Pix) {I-K') ' yW{y)Xj{y)) 



i,j>0 
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Applying this to Lp shows that 



Qk{x) = x^Q{x) - ^ {vjQ,{x)-UjP,{x)) 



i,j>0 



and applying it to ip shows that 



Pkix)=x''Pix)- J2 {wjQ^{x)-VjP^{x)) 



(2.12) 



(2.13) 



i+3=fc-l 



These are the recursion relations for Qfc, Pk- Taking the inner product of both sides of the first one with tpX , gives 



Vk=Vk- ^ (VjVt - UjWi) , 

i+3=fc-l 

ij>0 

recursion formulas which can be used to express the Vk in terms of the Wi, Vi, Wi 



Finally, using the definition of Uj in (2.5a), the fact 



7^X,,(2/)-(-l)'=%-a.), 



da, 



k 



and (2^) we find that 



dak 



= (-1)*^ {R{x,ak),x^ipXj) + {-lfQ{ak)a{Lp{ak) 
= (-l)*^' {p{x,ak),x'LpXj)qk 



= {-lfQj{ak)qk 
where q^k = Qj{ak) {qok = qk)- Similarly, 



dvj_ 

dak 
dvj_ 

dak 
dwj_ 

dak 



{-if Pj{ak)Q{ak) = {-ifpjkqk 
{-lf Qj{ak)P{ak) - {-ifqjkPk 
(-l)^Pj(afc)P(afc) = {-ifpjkPk 



(2.14) 



(2.15) 

(2.16) 
(2.17) 
(2.18) 



where pjk — Pj{ak) {pok = Pk)- From ( |2.12| ), (2.13) and (2.14) we recall that qjk and Pjk are expressible in terms of 
the qj , pj , Ui, Vi, and Wi . 



C. The Case m{x) = 1 



In this section we derive those partial differential equations that depend upon the differentiation formulas (1.4) in 



the special case 'm{x) = 1. We let D denote the differentiation operator with respect to the independent variable and 
recall that if the operator L has distributional kernel L{x,y) then 



(2J9) 



Using the differentiation formulas it follows that 



I- + 1-) K{x, v) = ^^""^ ^^y^ {^{xmy) + ^xMy)) 

dx ay / X — y 

B(x)-B{y) , C(x)-C{y) , , , , 
' ^ ^ ^^!^(a;)V'(y) + ^-^ifiixMy). 



x-y 



x-y 



(2.20) 
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Let us write 



(2.21) 



Then 



A{x) - A{y) 



x-y 



= ^ a^+k+ix^y^, etc. 



j,fc>0 



and we obtain the identity (recall (2_^) and ( 2.19| )) 

[D,K]= aj+k+i{x'^{x)v^i^{y) + x^^:{x)y^^{y))Xj{y)+ ^ p^+k+ix' ^{x)y''i:{y)Xj{y) 

j,k>Q j,k>0 

+ ^3+k+ix''p{x)y''^{y)Xj{y) - Y{-l)''K{x,ak)S{y - at) 

j,k>0 k 

from which it follows that 

(I - K)-'] = J2 ^j+k+i {Qjix) (/ - K')-' y''^iy)Xjiy) + Pj{x) (/ - K^^ yV(2/)x/y)) 
+ Y.p,+k+iPAx) {I ~ K'y\''^{y)Xj{y) 
+ l3+k+iQAx) {I ~Ky^ y%{y)Xj{y) 
- ^k)p{ak, y). 

We now use this last commutator to compute Q'{x) and P'{x): 

Q'{x)=D{I-K)-^^{x) 

^ {I - KY^ Dip{x) +\d,{I - Ky'^ Lp{x) 

= (/ - Dip{x) + ^ aj+fe+i {vkQ]{x) + UkPj{x)) 

+ Yf^j+k+i'^kPj{x) + ^7j+fe+i"fc<3i(a;) 

and similarly 

P'{x) = {I - Ky^ Dyj{x) + (wkQA^) + ^kPj{x)) 

+ Y Pj+k+iWkPjix) + Y lj+k+iVkQ]{x) 
^Y{-lfR{x,ak)pk. 
Finally we use the differentiation formulas (^]J) and representations ( 2.21 ) to deduce 

(/ - A^-' D^{x) = (^M^) + PiPii^)) ' 
(/ - K)-' D^{x) = - c^jPji^)) 

and so substituting into the above gives 



(2.22) 



Q'ix) = ^ + X! (^j+k+iVk + Y 7j+fc+lWfc ) Q]{x) 

j>0 \ k>0 



k>0 



'Y \f^^ OLj+k+lUk + Y Pj + k+lVk ) Pj{x) 

j>0 \ k>0 



k>0 



2m 



-^(-l)'=i?(x,afc)<7fc, 
fc=i 

j">0 y fc>0 k>0 J 

j>0 \ fe>0 
2m 



k>0 



k=l 



From (2.10), (2.11) and these last identities we deduce the equations 



dqi 



j>0 



k>0 



k>0 



'Y \^1~^Y ^J+k+l^k + Y (^J+k+l'"k I P]i 
j>0 \ k>0 k>0 

^(-l)'''i?(a„afc)gfc, 

k^i 



dp, 



j>0 



k>Q 



fe>0 



' Y + Y (^3+k+lVk + Y P]+k+lWk I Pji 

i>0 \ k>Q k>0 

^(-l)'=i?(a„afe)pfc. 

kyii 



Using (1^, (|2.23| ), and (|2.24| ) we deduce 



R{ai,ai) = X! +^aj+fc+it;fe + ^7j+fc+iWfc gj^pj 



fe>0 



fc>0 



j>0 \ k>0 k>0 

'Y M-J' ^ X! ^j+k+iU!k - Y 7j+fc+i«fc I 9j« 9i 

j>0 \ k>0 k>0 

i>o y fc>o fc>o 
^k{(l^Pk -PtQkf 



k^i 



ai - ak 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



We end this section with two differentiation formulas for R{ai,ai). From ( ^ ) and ( 2.22| ) we deduce that for x,y G J 
d d 5 \ 

^ + ^ + Ri^-Jj) = E«j+fc+i (Qj(x)Pfc(y) +Pj(a^)Qfc(y)) 
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' Pj+k+iP] {x)Pk {y) + 7j+fc+i Qj {x)Qk {y) - i-^f^i^^ ak)R{ak , y). 



Hence, using the chain rule, 



d 



—R{ai,ai)^2 ^ a^+k+i qjiPki + ^ Pj+k+iPjiPki 



j,k>0 j,k>0 



j,k>0 k^i 

A variant of this foUows from it and ( ^.9[ ) by the chain rule: 

—Rtj{tai,tai)^2ai ^ a^+k+i QjiPki + at ^ Pj+k+iPjiPki 
j,k>0 j,k>a 

+ai Y Tj+fc+igjj^fci +^(-l)''(afc - ai)-Rtj(^ai,iafc)" 

j,k>0 k^i 

Here the subscripts tJ indicate that this is the underlying interval. 



(2.28) 



(2.29) 



D. The Case of Polynomial m 

Now let us see how the above derivation has to be modified if m{x) is an arbitrary polynomial. In this section M 
denotes multiplication by m{x) and D continues to denote differentiation with respect to the independent variable. 
In place of the commutator [D, K] we consider the commutator 

[MD,K] = (m{x)-^+m{y)^+m\y)^ K{x,y)Xj{y) 



dx 



ak) 



while using ( |1.4D we compute that 



d_ 

' dx 



d m{x) — m{y) 
' dy 



X — y 



m{x)-^ + ■m{y)^ H — ^— ) K{x,y) = the right hand side of (2.20). 



(2.30) 



Therefore if m{x) is linear and we replace D by MD on the left side of ( 2.22 ), then the right side h as to be changed 
only by insertion of factors m{ak) in the last summand. It follows from this that ( 2.25| ) and (2.26) require only the 
following modifications: 



Insert on the left sides of (2.25) and (2.26) the factor m{ai). 
Insert in the last summands on the right sides of ( ^.25 ) and ( 2.26 ) 
the factors m{ak). 



while ( 2.27 ) and ( 2.28 ) require the following: 

Insert on the left sides of (2.27) and ( 2.28| ) the factor m{ai) in front of R{ai, ai) 



For general m(a;), if 



then 



Insert in the last summands of (2.27) and (2.28) the factors mlak) 

m{x) = Y i^kx'' 
m{x) ~ m{y)\ _ ^ ^ ^ 



(2.31) 



(2.32) 



1 



x-y 



m'[y) - 



x-y 



fij+k+2 x^y^ 



j,k>0 

j+k<dGg m—2 
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It follows from this and (2.30) that 



m{x)— — h m{y)— — (- m'{y) J K{x,y) = the right hand side of ( 2.2C| ) 



- ^(fc + l)fij+k+2 x'y^ {if{x)ip{y) - ip{x)ip{y)) 



So for m{x) of degree greater than 1 the right sides of ( |2.23| ) and ( ^.24 ) must also be modified by the addition, 
respectively, of the terms 

- '^{k + l)A^J + fe+2 (vkQjix) - UkPj{x)) , 

- + l)Aij+fe+2 (wkQjix) - VkPj{x)) . 

The upshot is that in this general case (|2|2|) and ^I^ ) require, in addition to ( |2.31| ), the following modifications 

Add to the right side of (2.25) 



j,k>0 

j + fc<dcg m— 2 



Add to the right side of ( |2.26| ) 
{k + l)^j+fc+2 {wkQjt - VkPji) ■ 



j,k>0 

j-\-k<dcg m— 2 



And for general m{x) we must modify ( 2.27 ) and ( ^.28 ), in addition to ( 2.32| ), by the following: 

Add to the right side of (2.27) 

- X! l)Mj+fc+2 {vkQji - UkPji) Pi 



j+k<~dcg m — 2 



j-\-k<.deg m — 2 



Add to the right side of (2.28) ^ (j — k)fij-^k+2Pji Qki- 



j,k>0 

j+A;<dcg m — 2 



(2.33) 



(2.34) 



(2.35) 



The identity ( 2.29| ) also requires modification if we do not have m{x) = 1, but since we shall only use it in this special 
case there is no need to write down the modification. 



E. The Exponential Variant 

Here we consider kernels of the form 

where b can be an arbitrary complex number. Because of the different denominator it turns out that the differentiation 
formulas should now be of the form 



(2.37) 



m{x)ip' {x) — (^A{x) + -'m{x)^ fi^) + B(x)ip{x), 
m{x)ip' {x) — —C{x)ip{x) — ^A{x) — —m{x)^ ip{x) 



where A{x), B{x), C{x) and 'm{x) are "exponential polynomials," finite linear combinations of the exponentials 
gfcbx ^-^ ^ _|_2^ . . .y -y^g compute, as the analogue of ( ^.30[ ), that 
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d 



Now, of course, we write 



^ ' dx ^^'dy 2 e^^ - e''y 

= the right hand side of ( 2.20| ) with denominator 
replaced by e^^ — e^. 



d b {m{x) - m{y)){e'"' + e'>y) 



K{x,y) 
x-y 



(2.38) 



A{x) 



etc.. 



the summations summing over negative and nonnegative indices, and 



A{x) ~ A{y) 



^jhx ^kby 

-k+ie-' e " 



E 

j,k<- 



aj+k+ie^'^e^'y, etc. 



(2.39) 



What arise now are functions Qk,- ■ ■ , Wk defi ned, for ne gativ e as well as nonnegative values of k, by replacing x^ by 



e in the earlier definitions. Analogues of ( 2.12 ) and ( 2.13 ) hold for negative as well as positive values of k so all 



the Qk and are expressible in terms of Qo, Po, as well as the Vk in terms of th e u^., w^, Wk- 

Notice that if 771(2;) is constant t hen th e third term in the large parentheses in (2.3S) vanishes and so we obtain in 
the end the analogue of ( 2.25 ) and ( 2.26 ); in addition to the change in the range of indices now and the fact th at the 
double sums have two parts, as i n (2;3S| ), we must in the single sums over j add to the terms aj and P-j in ( | 2.25| ) 
and the terms —jj and —aj in ( 2.26 ). For general m{x) we must insert factors m{ai) on the left sides of ( 2.25]) and 



(2.26) and factors m{ak) in the last summands on the right, and then add terms coming from the difference 



■m'{y) 



h {m{x) - m{y)){e'"' ^ e^y) 



'_(p{x)'il;{y) - 'ip{x)ip{y) ' 



(2.40) 



as at the end of the preceding section. We shall not write these down since in the only case we cons ider later we have 
m{x) — 1. Two of the equations involvin g R{ x, y) must also be modified. We see first that in ( 1.7a| ) the denominator 



must be replaced by e 



^bak 



Second, 



must have the factor h e inserted on the left side, with the result that 



( ^.27 ) must have the factor foe'"'' inserted on the left side. Note that ( 2.29 ) is unchanged 



Remark 1 . The product of the first two factors in (2.40 ) is a n exponential polynomial in e^^ and e^y . It was precisely 
to achieve this outcome that we required the formulas ( ^.37 ) to have the form they do. 

Remark 2. In case 6 is real the change of variab le x ^ transforms the operator with kernel (|1.2|) acti ng o n the 



set e to an operator with kernel of the form (2.36) acting on J, with the new (</?,?/') P&ir satisfying (2.37). So 



we see that there is m ore than simply an anology between the two situations. In fact we could have allowed the 
various coefficients in ( |1.4| ) to be linear combinations of negative or nonnegative integral powers of x, and then the 
two situations would have been completely equivalent for real h. 



III. SINE, AIRY AND BESSEL 
A. Sine Kernel 



The simplest example is the sine kernel 



where we take 



K{x,y) 



A sin(a; — y) 
TT X — y 



ip{x) = \j — sinx, ij^ix) — y — cosx. 

The differentiation formulas hold with 

m{x) = 1, A{x) = 0, B{x) = 1, C{x) = 1. 
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( It i s useful to incorporate a parameter A G [0, 1] into K] c.f. for mula (1.1).) The partial different ial eq uatio ns are 
(1.8), ( p^ ) (the universal equations along with universal relation ( 1.7a )), and the specialization of ( 2.25 ) and ( 2.26| ) 
which now read, respectively. 



dai 

dpi 
dai 



Pt - '^{-l)''Ria„ak)qk , 



-qi - ^i-'^)'' R{ai, akjPk , 



(3.1) 
(3.2) 



along with the specialization of ( 2.27 ) 



i?(a„a,) =p,2 + (zf + ^(-1) 

k^i 



kili Pk -Piqk) 

CLi - ttk 



(3.3) 



These are the equations of JMMS |^ though they appear here in a slightly different form due to our use of sines 
and cosines in the definitions of (p and -0 rather than the alternative choice of e^*^, which we could have taken just 
as well. (They also appear slightly different in p7| due to our convention here not to put in a factor of tt into the 
definition of the sine kernel.) 

For the case of a single interval J — {~t,t), s — 2t, these equations imply that <t{s;X) := —sR{t,t) satisfies the 
Jimbo-Miwa-Okamoto a form of Painlcvc V. We refer the reader to the literature for a derivation of this, properties 
of the solution of this equation, and the implications for random matrices ||,|l^,^,|9|,|^,|7|j 



B. Airy Kernel 

For the Airy kernel we have (again inserting a parameter A into K) 

ip{x) = VAAi(x), Tpix) = \^Ai{x) 

from which it follows that 

m{x) ^ 1, A{x) ^ 0, B{x) = 1, C{x) = -x 

since Ai"[x) = a;Ai(a:). For notational convenience we write u = uq and v — vp . In a ddition to the universal relations 
( 1.7a )-(1.9), we have two additional equations for u and v, viz. ( 2.15 ) and ( 2.16 ) for j — 0. Using the recursion 
relation ( 2.12| ) for k = I, we deduce that (2.25) and ( ^.26 ) reduce to 



-uqi 



Pi - ^[-lfR{a^,ak)qk , 



k=^i 



dqi 
dai 

(a^ - 2v)qi + up., - ^{-l)''R{a,, ak)pk 
<ja,i . 

k^i 



(3.4) 
(3.5) 



and (2.27) reduces to (again using the recursion relation ( p.l2| ) for fc = 1) 

D/ \ 2 2 f, , o 2 , V^/ ^-.kiliPk ^ PiqkY 

R[ai,ai) ^Pi ~ aiqi - 2uqip, + 2vqi + } (-1) 

^ at - ak 



k^i 



(3.6) 



These are the equations derived in | |3^ . We mention that in addition to these equations, two first integrals were 
derived which can be used to represent u and v directly in terms of the qj and Pj (see (2.18) and (2.19) in Q). We 
also remark that in the case J — (s,oo), the quantity R(s, oo) was shown to satisfy the second order nonlinear a DE 
associated to Painleve II. Again we refer the reader to [B8| for details. 
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C. Bessel Kernel 



For the Bessel kernel 

from which it follows (using Bessel's equation) that 

m{x) = X, A{x) = 0, B{x) = 1, C{x) = ^{x~a^) . 
Again using the recursion relation ( ^.12 ), we deduce that (2.25) and ( 2.26| ) become, with the additional insertions 



UiT^ = -UQi + Pi - ^[-l) akR{ai,ak)qk 



da, 4' 
dp 1 1 



and ( 2.27 ) with insertions ( 2.32 ) becomes 



a,Ria„a,) = {a^ - a, + 2v) q, + ^uq.p, + + V(-l)'=afc^^^^^-^^ . 

(As before, u — uq and v — vq.) These are the equations derived in ||3^. As was the case for the Airy kernel, two first 
integrals were derived which can be used to express u and v directly in terms of the g, and pj. For the case J — (0, s), 
the quantity cr(s) — si?(0, s) was shown fs^] to satisfy the a DE for Painleve 111 [p4 35|. 



IV. BEYOND AIRY 



In this section we give as an example of our general system of partial differential equations the simplest case "beyond 
Airy" in the sense discussed in the Introduction. In the language of 2D quantum gravity matrix models (see Q for a 
review), we are considering the case of pure gravity. Thus we take 



= Dl + -ax)D^ + -i'{x) 

where Dx is differentiation with respect to x. The string equation implies that ^(x) satisfies 

e"'(a;) + 6f(x)^'(2;) + 4 = 
which when integrated is Painleve 1 [^ UT^j2l[| : 

i"{x) + U''{x)+^x = 0. 

(Without loss of generality we may set the constant of integration to zero since it corresponds simply to a shift in 
the variable x. And, of course, the '3' and '4' can be changed by scale transformations to give the canonical form of 
Painleve I.) Exactly what solution S^{x) one choose s for pure gravity is still of some debate (on this point see ill|| an d 
references therein). The function ip{\,x) satisfies ( |l.ll|) and ( 1.12 ) which implies that if we define il){\,x) by ( |l.l3[) , 
then the differentiation formulas are 

m(A) = l, A{\)^~\£,'{x), B(A) = A+ic(^), 

c{\) = + ^ax)x + h e{x) + i r 
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where we remind the reader of the change of notation in the independent variable (see Introduction). 

Sin ce C(A ) i s qua dratic in A, the equ ations will involve m^, Vj, and Wj for j = 0, 1. Using the recusion relations 
( ^.12D , ( |2.13| ), ( |2.14D the equations ( |2.25| ) and { \2.26\ ) specialize to 



d 



-^(-l)"-i?(a„afe) 



9A; 



where M{ai) is the 2x2 matrix whose elements are given by 

Mii{ai) ~ ~ w + -£,u + uv — ui — a.iU, 

7Wi2(a,) = ai-u'^ + 2v + i^, 

= - - I ^" - ^Oi-y + - 2wi -2uw-^ a,^ + 
Al22(a») = -Miiia,). 
Similarly, ( 2.27 ) specializes to 

I ^ n ■ — n I. 



The universal equations are (|1.7a))-(p^) and (|2.15|)-(^l8|). 

For the case J — {s,oo) (this should be compared with the analagous case in Airy |38)), we are able to find two 
first integrals that allow us to eliminate the quantities ui and vi. (It is natural to take the boundary condition that 
all quantities evaluated at 00 vanish.) We denote hy q = q{s), etc. the quantities corresponding to the first endpoint 
fli = s. The first relation is quite simple 

q'^ ~ 2uui + + 2vi + i^-it^ - ^£,'u = 0, 

but the second one we found is rather messy 

— + sq^ + u + 2pqu + xu^ + q^u^ — 2u'^Ui — u\ — 2xv 
—Aq^v + Quuiv + "iu^v^ — St;"^ — 2u'^w — 2uxW + ^uvw — vo^ 
1 



-L 13 

-q^ + —u^ + 2uui — 6v^ + 2uw — ^'^^^^ ^ 
/ 1 , 1 3 1 1 







where q, p, u, etc. have argument s and ^ has argument x. (The variable x appears since we used Painleve I to 
eliminate ^"(x) in our equations.) 

Using these first integrals to eliminate ui and vi (note that w also drops out) we obtain the system of equations: 



1 



1 



8- 
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xs + s'^ + -x^ + -^-^ + —(f) ] q+p+ (As- iv- C)q + "^uq P - 2gp^ 



1 



1 



1 



p" = [2s^2v- -n q + [xs + s-" + -xi + + — (f) +2u]p 
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+ (4s - 4w - ^) pq^ + Auqp^ - 2p^ 
and, of course, we still have the universal equations 



-qp. 



Letting R{s) = R(s,oo) we find from ( |2.28| ) and ( |2.35| ) that 

1 



R' 



2s + 2v+^i]q^ +p^ - 2uqp. 



Clearly, further analysis of these equations is needed to be able to analyze the associated Fredholm determinant. For 
example, can one derive a differential equation for R itself? 
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V. FINITE N HERMITE, LAGUERRE, JACOBI AND CIRCULAR 



A. Hermite Kernel 



1. The partial differential equations and a first integral 

It follows from the Christoffe l-Da rboux formula that the kernel XKN(x,y) for the finite N Gaussian Unitary 
Ensemble (GUE) is of the form (L2) provided we choose 



with ipk{x) the harmonic oscillator wave functions 

(^fc(x) = -=L=e--'/2iJ,(x), fc = 0,l,... 

and Hk{x) are the Hermite polynomials iQ. This is well-known and we refer the reader to for details. From the 
differentiation and recurrence formulas for Hermite polynomials it follows that the differentiation formulas for cp and 
ijj hold with 

m(x) = 1, A{x) = -X, B{x) = C{x) = \/2N . 
We can therefore immediately write down the equations 



1/4 



daj 
dpj 



a,p, ~ (Vm + 2w) q, 



ttj — ttk 



Pk 



along with 



R{a,,a.^ = -2a,p,q, + (V2N - 2u) p] + (Vm + 2w) g| + (^Z^L_M^il 



(5.1) 
(5.2) 

(5.3) 



d 

—R{aj,aj) = -2pj qj - '^{-l)''R{aj,akf , 



(5.4) 



^Rtj{taj,taj) = -2ajPjqj +^(-l)''(afe - aj)Rt j{taj,takY 



These follow from formulas (p.25|)-(|2.29D. 

We now derive a first integral involving u, w, pj and Qj] namely we show 

V2N {u-w)+ 2uw = - ^{-lypjqj. 



(5.5) 



(5.6) 



Observe that 



(E a^) P^l^ = + 2-) (V2iV - 2u) 



and 
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Multiplying the first equation by (—1)-' and summing both equations over j, results in 



It follows that the two sides of ( p.6[ ) differ by a function of (ai, . . . , a2m) which is invariant under translation by any 
vector (s, • • • , s). Since, clearly, both sides tend to zero as all — > oo, their difference must be identically zero. 



2. Bulk scaling limit of finite N equations 



We now show how (p.l|)-( |5.3| ) reduce to the sine kernel equations ( |3.1[ )-(p.3|) in the "bulk scaling limit." For a fixed 
point z, i.e. independent of iV, the density p{z) in the GUE is asymptotic to \/2N /-k as N ^ oo. The bulk scaling 
limit corresponds to measuring fluctuations about this fixed point z on a stretched length scale proportional to \/2N 
and then taking N — > c». Denoting for the moment the bulk quantities with a superscript B, this means we set 




and consider the limit — > oo, — > z such that is fixed and 0(1). In this limit we deduce from the asymptotics 
of the harmonic oscillator wave functions (see, e.g.. Appendix 10 in |28|) that both and ij) are 0(1) quantities in the 
bulk scaling limit. From this and the fact that it is K{x, y)dy which is 0(1), we deduce that both qj and pj are 0(1) 
quantities in the bulk scaling limit. An examination of the inner products defining both u and v shows that these too 
are 0(1) quantities. Thus if we formally replace 



z + 



9i 



Pj 



R{aj,aj) ^ V2N R''{af,af) 



in (5.1)-(5.3) (and replace all derivatives by derivatives with respect to af), take ^ oo, we obtain (p.l|)-(3.3) 



3. Semi-infinite interval and Painleve IV 



In this section we specialize the finite N GUE equations to the case of m = 1, oi = s and 02 = 00, i.e. J = (s, 00). 
We write q{s), p{s), and R{s) for gi, pi, and i?(ai, ai), respectively, of the previous section. The differential equations 
reduce to (' = d/ds) 



q = -sq 
p' — sp — 



V2N - 2uj p, 
2N + 2w^ q, 



u 
w' 



-p 



(5.3) reduces to 



and the first integral is now 



i?(s) = -2spq + V2N - 2ii + V2N + 2w]q 



'2N (u — w) + 2uw = pq. 



(5.7) 

(5.8) 

(5.9) 
(5.10) 



(5.11) 



(5.12) 



We proceed to derive a second order differential equation for R{s) and show that it is a special case of the Jimbo- 
Miwa-Okamoto a form of Painleve IV |24|j3^]. Relation (5.4) is now 
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R' = -2pq, 



(5.13) 



while (5/7) and (5^) give 

{pqY = (Vm - 2uj - (Vm + 2u?j . 

Differentiating one more time gives 

(pq)" = 2s - 2uJ p^ + (^V2N + 2w^ q^ - 8pq (^N + V2N {w - u) - 2uu?j + ip^q^ 

I (Vm - 2u^ p^ + (Vm + 2w^ g2 1 _ gjsipq ^ i2p2q2 



^ 2s- 



where we used the first integral ( 5.12 ) to obtain the seco nd eq uality. Referring back to ( 5.11 ) we see that the term in 
curly brackets in the last expression is i? + 2spq. Using ( ^.13 ) to eliminate all terms involving pq in the last equation 
we find 

R'" = -4s{R - sR') - 8NR' - Q{R'f. 
This third order equation can be integrated (the constant of integration is zero) to give 



{R"f + A{R'f {R' + 2iV) - 4 {sR' - Rf = 0. 



(5.14) 



Comparing this with (C.37) of [g^ (see also |Q]), we see immediately that this is the a version of Painleve IV with 
parameters (in notation of |2^]) vi — Q and V2 = 2N . Explicitly in terms of the Painleve IV transcendent y = y{s) 
we have 



(5.15) 



with Piv parameters a = 2N — 1 and (i — Q. Recall that w = wlz) is a Painleve IV transcendent with parameters a 
and /3 if it satisfies the P/y equation 



1^ 



1 



2w V dz 



dw 



3 5 

+ -w^ + Azu? + 2(z^ - a)w + — 
2 w 



We are, of course, interested in the family of solutions that vanish as s — > +oo. 

This particular Pjv has been studied by Bassom, et. al. ||] (see also |^^). To make contact with their notation 
define 



{y is the above Piv transcendent) so that rj satisfies 



3rf + 2^77^ 
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with ly = N — I. They analyze the one-parameter family of solutions rjkiS,]!^) satisfying the boundary condition 
ry(cx)) = 0. The parameter k is defined uniquely by the asymptotic condition: 



?7fe(C; ~ fc^^exp 



In terms of our parameter A we have 



4 



A 



as C 



oo. 



23/2(iV-l)!0f 



(This identity is derived by examining the large positive s asymptotics of i?(s), which is easy because of the rapid 
decrease of the kernel as s ^ +00.) These authors prove that for all positive integers TV > 1 the solution 77fe(^, — 1) 
exists for all ^ whenever A < 1, and that rik{(,,N — 1) blows up for a finite ^ whenever A > 1. These results are in 
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complete agreement with what one expects from the spectral theory of the Fredholm determinant. There are formal, 
but not rigorous, results that solve the connection problem for the asymptotics as ^ ^ — oo; in particular, for A = 1 



1/2 



as ^ —oo. 



(5.16) 



Using (5.16) (and computing higher order terms by using the differential equation) we find that 



R{s) = -2Ns + ^ - 



N iV2(l + 9iV2) iv3(io + 277V2) 



4s5 



+ 



4s7 



+ 



as s 



(5.17) 



4-. Distribution functions for X,nax and Xmin 
If we denote the smallest and largest eigenvalues of a matrix from the GUE by Amin and Amax, respectively, then 



in the notation of (1.1) we have 

P(Amax < s) = E{0; (s, oo)) = det (/ - K). 
Thus using (|l.3| ) we deduce the representation 

P(Amax < s) = exp | - ^ R{t-, 1) dtj 

where i?(s; A) denotes the function R{s) of the preceding section with parameter value A. This is our representation 
of the distribution function for Amax in terms of a Painleve transcendent. There is of course a similar representation 
for the distribution function for Amin- 

The authors of PpO| give an algorithm to compute the quantities rikiS,, v), v — positive integer, of the last section 
exactly in terms of the error function 

X2 



/(C) = 2 exp(^-yj dx. 

That such elementary solutions of the P/y transcendent exist, at least for the case A = 1, is now clear from the 
random matrices point of view since E{Q\ (s, oo)) is expressible in terms of integrals of the form 



This follows from (1.15) with / the characteristic function of (— oo, s). 



5. Edge scaling limit from Painleve IV equation 
The edge scaling limit |38| corresponds to the replacements 

s ^ V2N + , and R V2N^'^R 

V2iVi/6 



m dsl^ ) and retaining only the leading order term as — ^ oo. The result of doing this is 

{R"f + 4 {R' f + 4i?' (i? - sR') = (5.18) 



which is the equation derived in ||3^. We remark that (5.18) is the cr form for Painleve II, see (C.17) in ||2J] and 
Proposition 1.1 in |M. 
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6. Symmetric single interval case 



In this section we specialize the finite N GUE equations to the case of m = 1, ai = — t and 02 — t, i.e. J = {—t, t). We 
denote by q{t) and p(t) the quantities (72 andp2, respectively. Since (^Ar(— — {~l)^ifi^{x) and K{—x, —y) = K(x,y), 
we have qi = (—1)^(7 and pi = — (— We further set 



and record that 



R{t) R{t, t) = -i), R{t) := (-l)^i?(-t, t) = {-lfR{t, -t) 

dlogD{J;X) 



dt 



-2R{t). 



Now (p is even or odd depending on whether N is even o r odd , with having the opposite parity. It follows from this 
fact, and our choice of sign in the definition of i?, that ( 1.7a ) specializes in this case to 



R{t) 



' t 



while (5.3) and (5.5) specialize to 

R{t) = -2tpq + (Vm - 2u) / + (Vm + 2w^ q^ - jp^q^ 



(5.19) 



(5.20) 



^ = 2tR + 2R' 
dt 



(5.21) 

The last is the finite N analogue of the Gaudin relation. (See, e.g. [E8^,p7|.) The differential equations specialize to 



dq 
It 
dp 
'dt 

du 

Hi 
dw 



2Rq-tq + (V2N -2u]p, 



dq2 ^ dq2 
dai da2 

-2Rp + tp- ( V2N + 2w]q, 
2p'. 



(5.22) 
(5.23) 
(5.24) 
(5.25) 



And the first integral (5.6) is now 



/2iV (u-w) + 2uw = -2pq. 



(5.26) 



7. Differential equations for R and R 



It follows easily from ( |5.19| ), ( |5.20D and (|5.2l| ) that 

d 



dt 
d 

di 



{tR) = ^{r'), 



At^R 



where 



'2N + 2w + ip^J V2N - 2u 



(5.27) 
(5.28) 



Equations ( |3.27 ) and (5.28) are the finite N analogue of those derived by Mehta (See also discussion in ^7\.) 

We now eliminate the quantity r. For this derivation only, we write a(t) := tR{t) and b{t) := tR{t). We begin with 
the obvious 
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|r|4 = 5R(r2)2 + 3(^2)2. 



(5.29) 



Now 



3(r2)2 = 4p2g2 ^^/2N - 2uJ (Vm + 2w 

= Ap^q^ {2N + Apq) 
= 4^2 {2N - 4h) . 



(5.30) 



We now use ( ^.27 ) and ( ^.2S| ) to obtain expre ssions for 5R(r2)2 and |r|'^, respectively, and the above identity for 
5(r2)2. These expressions, when used in (5.2£), we have an equation for a, 6, and their first derivatives. If we use 



the generalized Gaudin relation (5.21) to eliminate the appearance of da/dt (the one appearing to the first power), 
we obtain 



8a6 + 8iV&2 



and together with (5.21), which in the a and h variables reads 

' da^ 



t 



dt 



a + 2^2 + 2i2&, 



(5.31) 



(5.32) 



we have two differential equations for R and R. 

Eliminating a, we obtain a single second order equation for b and therefore R: 

[tR!' + 2R' -24:t^ R^ + 8NtRy ^ 4: (2R - ^ (^St'^R^ (^N - 2tRj + {r + tR'^ ^) = 0- (5-33) 

This last equation is the finite N analogue of (1-18) of Mahoux and Mehta We could, in a similar way, eliminate 
b and so obtain a second order equation for R, but the result is messy and we shall not write it down. 



8. Small t expansions of R and R 



The boundary conditions at i = for ( |3.2l| ) and (5.33) follow from an examination of the Neumann expansion of 
the resolvent kernel. Setting pq :— K{0, 0), the density of eigenvalues at 0, we find 



and 



Rit) = po + 2plt + p„ (4p2 + (-i)A^) t' + 2^ (9p2 _ 27V + (-1)^') t3 + 

^ {288p^o + 40p2 ((-1)^ - 2A^) - 12(-1)^A^ -3)t^ + 0{f) 
18 



i-irm) =Po + 2plt + ^ (I2p2 -4N- i-ir) t' + ^ {9pl -2N+ i-lf) t 



N\ .3 



Po 
90 



(I440p^ + 200p^ ((-1)^ - 27V) + 487V2 + 12(-l)^iV + 9)t'^ + O(t^) 



9. Level spacing probability density pN{t) 



(5.34) 



(5.35) 



For m = 1 if we let EN{0;ai,a2) denote the probability that no eigenvalues lie in the interval (01,02) and 
p^f (0; fli, 02) da2 the conditional probability that given an eigenvalue at ai the next one lies between 02 and 02 + da2, 
then the two quantities are related by 



PAr(0; 01,02) = 



1 92^^(0-0,^ ^2) 
p(ai) daida2 
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where p{ai) is the density of eigenvalues at ai. From the expression for the logarithmic derivative of the determinant 
(with A = 1), we have 



9^£;Ar(0;ai,a2) _ dR{ai,ai) 



daida2 



da2 



En - R{ai,ai)R{a2,a2)E^ 



N 



Differentiating (5.3) (with m = I and j = I) with respect to 02 we obtain 

2 



which when evaluated at ai 



dR{ai, ai) 
da2 

~t, a2 — t becomes 

dR{ai,ai) 



qiP2 - Pig2 
ai - 02 



da2 



~\ 2 



ai — — t,a2—t 



Calhng pN{t) ■— Pn{0', —t, t) we thus obtain 



Pit) 



-R\t)\EN{t) 



where E^it) = En{0; —t, t) (we used also p{—t) = p{t)). 
Using the expansions of i?, R and En [t) we find 



PN{t) = 3 (A^ + (-1)'') P^t" - ^ (16A^' + 2%-lY'N + 13) Pot" + 

— (1287V^ + 452(-l)^iV2 + 5297V + (-1)^205) pot^ + ■■■ . 
315 ^ ' 

Not only does this hold for fixed N and but it also holds uniformly in N and t as long as t = 0{N^^/'^). The reason 
is that in this range of the parameters the operator K has norm less than a constant which is less than 1 and has 
bounded Hilbert-Schmidt norm. Thus the Neumann series for the resolvent kernel converges in trace norm. 
To compare with the bulk scaling limit we replace t by t/poi and deduce (recall po ^ V2A^ /tt as N ^ 00) that 



Pit) 



1 



t 



2^4 



N^oo pq pq 3 45 315 



which is the well-known result Observe that the large N corrections to these limiting coefficients are 0(1/A^). 
(Note that we inserted a factor of n in our definition of the new t variable so as to have the same normalization as in 



B. Laguerre Kernel 

1. The partial differential equations 

Again by the Christoffel-Darbo ux fo rmula it follows that the kernel for the finite N Laguerre Ensemble of N x N 
hermitian matrices is of the form (1.2) provided we choose 

(p(x) = VA (iV(iV + a))'/ViV-i(^), ijix) = V\iN{N + a))^^\Nix) (5.36) 
(cat in (^^) is negative now) where 



and L'^{x) are the (generalized) Laguerre polynomials See Chap. 19 of and for further details and refer- 
ences. Fr om t he differentiation and recurrence formulas for Laguerre polynomials it follows that we have differentiation 
formulas (|l.4|) for ip and ip with 



23 



I n 

m{x)^x, A{x)^-x---N, B{x) = C{x) = y/N{N + a) . 



We therefore have the equations 



1 a 

2 ^ 2 



Uj R{aj,aj) = (flj -a-2N) q^pj + (^^yN{N + a) + p"^^ + (^^N{N + a) - wj q: 



flj — ak 



-aj R{aj,aj) = qj pj - ^(-l)''a/c-R(aj, Ofc)^ 



These foUow from formulas (^)- (|2!28|) as modified by (|!3|) and ( ^!32| ) 



(5.37) 
(5.38) 

(5.39) 
(5.40) 



2. Single Interval Cases (0, s) and (s, oo) 

We c onsid er first the interval (0, s). We set ai = 0, 02 = s,q2 — q,P2 — P,R{s,s) = R{s) and find that eqs. 
( ^.37D ~( |5!40| ) with j = 2 speciahze to 



sq' ^ Qs - I - Tvj 9 + (^^N{N + a)+ujp, 

sp' = -(^^N{N + a)-w'^q- (^s-f 
s R{s) (s - a - 2N) qp + (^y^N{N + a) + uj p^ 

+ (V^(iV + a) - w) q"^, 
{sR{s)y = qp 



(5.41) 

(5.42) 

(5.43) 

(5.44) 
(5.45) 



(notice that the t erms corresponding to fc = 1 in the double sums on the right sides of (5.37)-(5.40) are equal to zero) 
while ( 2.15 ) and ( 2.18 ) specialize to 



u' = q'^, w' = p^ . 



Tedious but straightforward computation using (5.41)-(5.46) gives 

s{pqy ^ [^/N{N + a) + u) p^ - [^/N{N + a) - q^, 

s^ipq)" = (2iV + a - s) { (^,/N{N + a) + p^ + (^^/N{N + a) - u-) q^ } 
-|(V7V(iV + a)+M) (^^Ar(iV + a) - g^j +2sp^q'^ 
-AN {N + a)pq + ^uw + \/N{N + a) (w - m)| 4pg. 
Now it follows from ( ^.45| ) and ( ^.47| ) that 



(5.46) 



(5.47) 



(5.^ 
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s R{s) — sqp 



(5.49) 



has derivative 



But it follows from (5.46) that 



UW+ \/N{N + a) {w - u) 



has exactly the same derivative. Hence the two must differ by a constant. This constant must be since (5.49) clearly 
vanishes when s = 0, and so do u and w. Thus we have derived the identity 



t w + y/N{N + a) {w~u) ^ s R{s) 



sqp. 



(5.50) 



Now we ca n see that that every term in (5.48) can be expressed in terms of R{s) and its derivatives (up to order 
3). By (5.45| ) this is clear for all pr oduct s pq and its derivatives. This is true of the first e xpression in curly brackets 
in (|5.48|) by what wc just said and ( 5.44 ), of the first expression in curly brackets by ( 5.47 ), and of the last expression 
in curly brackets by ( ^.5C| ). 

Thus we have derived a third-order differential equation for i?(s). In terms of (t(s) :— s R{s) it reads 



s'^a'" 



(2N + a-s)(j- 



4:Ns-2as)a' -sa" + 6s {a'y ~4aa' 



It follows from this that the two sides of a purported identity 



(S(T' 



ll\2 



4 s {a'y + {2a + 4N-2s)aa' 
+ {a'^ -2as-4Ns + s^) {a'f -Act {a') 



(5.51) 



(5.52) 



differ by a constant. (The third-order equation is equivalent to the two sides' here having the same derivative.) Now 
it is clear that if a i s suffi ciently large then a is twice continuously differntiable up to s = and cr(0) = cr'(O) = 0. 
Hence both sides of ( 5.52 ) vanish at s = and so the difference in question equals 0. Thus the identity is established 
for a large. But both sides of the identity are (for s > 0) real-analytic for a > — 1 and so if they agree for large a 
they must agree for all a. 

Comparing ( 5.52| ) this with (C.45) of ||2^] we see that — (t(s) satisfies the a version of Painleve V with parameters 
Vf) ~ vi — 1^2 — N and u^^ — N + a. 

The boundary condition at s = for a depends, of course, on the parameter A in ( 5.36| ), and we write cr(s; A) 
instead of ct(s) to display this dependence. With the help of the Neumann series for the resolvent kernel, we compute 
the small s expansion 



ct(s;A) = Acos"+i [\ 



a 



2N 



a 



+A 



2 '--0 ^2a+2 



+A3 



1 - 

3a+3 



(2a + 3)(a + 2iV) 



(2 + a)2 



+ •••) 



where 



Co 



T{N + a + l) 



T{N)T{a + l)T{a + 2) 
And now, as in the case of the GUE, we have a representation 

^'(Amin > S) = exp |- / dt 



(5.53) 



There are only minor changes required in the above analysis when we take J = (s, oo) and this leads to an analogous 
representation for the distribution function of Amax- 
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3. Singular values of rectangular matriees 



If j4 is an X M rectangular matrix [N < M) whose entries are independent identically distributed complex 
Gaussian variables with mean and variance 1 then the N x N matrix A A* (whose eigenvalues are the squares of 
the singular values of A) belongs to the orthogonal polynomial ensemble associated with the weight function 

w{x) = x''-^ e~-/^ 

(See [0, Cor. 3.1.) It follows that the distribution of the smallest singular value of A is given by the right side of 
( p. 53 ) with s replaced by and, of course, a = M — N. There is a similar representation for the distribution 



function of the largest singular value of A. 



C. Jacobi Kernel 



The situation here is so similar to the preceding that we shall only indicate the main points. For the finite N Jacobi 
ensemble 



wix) = (l-a;)"(l + a;) 



/3 



with a, P > —1 and we must take in (1.2) 



ip{x) = y/XuN Vn{x), ip{x) = \/\aN ipN-l{x) 



where 



N{N + a){N + I3){N + a + 13) 



N+{a + l3)/2 y {2N + a + P- 1){2N + a + /5 + 1) ' 



From the differentiation and recurrence formulas for the Jacobi polynomals we have the differentiation formulas (1.4) 
with 



where 



ao 



ai 



70 



m{x) = 1 — x^, A{x) = ao + ai x, B{x) = (io, C{x) = 70, 



2{2N + a + (iY 
a + (3 



N - 



2N + a + (3 
2 

2N + a + l3 



y^N{N + a)iN + (3){N + a + 13) 



y/N{N + a){N + I3){N + a + l3) 



2N ^ 


- a - 


^13 + 


1 


2N ^ 


- a - 




1 


2N ^ 


- a - 




1 


2N ^ 


- a - 


h(3 + 


1 



Usi n g thes e the interested reader could without difficulty write down the Jacobi analogues of the general equations 
(^!37|)-(i5.40|) in the La guerre ensemble. We shall restrict ourselves he re to t he cas e J = (— 1, s), the analogue of the 
interval (0, s) in Laguerre, and find the following Jacobi analogues of (^)~(|5l6|): 



(1 — s^) q' = (ao + ais + v) q+ (/3o + {2ai — 1) u) p, 
(1 - s'^)p' = (-70 + {2ai + l)w) q- (ao + ais + v) p, 
(1 - s^)R{s) = 2 (ao + ais + v) pq + A + B, 
{{1 - s^) R{s))' -2aipq, 
v' = pq 



(5.54) 
(5.55) 
(5.56) 
(5.57) 
(5.58) 



where 
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A=iPo + (2ai - 1) u) S = (70 - (2ai + l)w) 



From ( |5.57D and (|5.58D we deduce 



(1 - s2)i?(s) ^2aiv 



(note that both sides here vanish when s = —I) and ( 5.54 ) and ( 5.55 ) give 

{l-s'){pqy = A-B. 

Another differentiation gives 

(l-s2)(pq)"^_2 (ao + ais + v) {A + B) + 2 s {A - B) 

+4ai (1 - s^)p^q^ - —AB. 

pq 



(5.59) 
(5.60) 



(5.61) 



Now (^.57 ) and ( ^.5£ ) show that pq and w are expressible in terms of R and its first derivative. And so, using (15^5^ ) 
and ( 5.6Cl| ), we deduce that A and B are expressible in terms of R and its first two derivati ves. F inally, using (5.61), 
we obtain a third order equation for R. Instead of this we write down the analogue of (|5.5l| ). If we define now 
a{s) := (1 — s^) R{s) then we have 



(1 - s^f a'" ^ {1 ~ 



2a\ 



2s(l-s")CT"-2(l-s2)(a')2 



-2(1 J- I - \ (ao + ai s) a. 



(5.62) 



For the boundary condition at s = — 1 we compute the small s + 1 expansion to be (with an obvious notation) 

, . , M+fl / /3(1 - a) - 27V(a + /3 + iV) , , \ 
a{s- A) = Aco(l + s)i+'^ [\ + '-^ 2{2^(5) + + ' ' ' j 

I A^cg ^2+2/3A I 2 + 6/3 + 3/3^-(3 + 2/3)(2jV^ + 2iVa + 2iV/3 + a/?) \ 

^2(l + /3)^ ^'^ 1, ^ 2(2 + /?)2 ^ '^-'^ J 



3^3 



^4(l + /3)2^^'*^ ^ + 



where 



Co 



r(A^ + g + /? + l)r(jV + /3 + 1) 

2^r(/3 + i)r(/3 + 2)r(iv)r(iv + a) ' 



We have not been able to find a first integral for (5.62), in other words a second order equation which is analogous to 
(P2|). 



D. The Circular Ensemble 

1. The partial differential equations 

If U (N) denotes the group of N x N unitary matrices, then the finite N circular ensemble of unitary matrices 
(sometimes denoted CUE) is this set U{N) together with the normalized Haar measure. Just as for the orthogonal 
polynomial hermitian matrix ensembles, the level spacing distributions arc expressed in terms of a Fredholm determi- 
nant of an integral operator defined now on the unit circle. All of this is well-known and we refer the reader to either 
Dyson or Mehta (2|] for details. 

The integral operator for the finite N CUE is 
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^ ^ ^ sinf (a; - y) 
27r sin i(x — y) 



where 



Thus the differentiation formulas ( 2.37 ) hold with b = i and 



m{x) = 1, ^(a;) = 2 ^' ^^^^ " '-^'^^^ " ^ 



and so from the considerations in Sec. II B we deduce that 



and 



R{aj,aj) ^iNpjQj - ^ {-if R{aj,ak) (qkPj -PkQj), 



dq^ N+l 



dpj _ .N~l 



J2{~l)'R{a„ 



a-kjPk- 



(5.63) 
(5.64) 



(5.65) 



2. Single interval case 



Now we specialize to J = (— i, t) and take oi — ~t, a2 — t. Write p, q for p2, 92 and note that since f{x) = ip{—x) 
and K{x, y) is even we have pi — p, qi — q- Let's also write 

R{t) R{t,t), R{t) R{~t,t) = R{t,~t). 



Then ( |2.29D with J = (-1, 1) and i = 2 gives 

R' = 2R^ 

O giA 



di 



nves 



whereas (5.63) with j — 2, k = I gives 



R{t)^Npq-lR{t) (qp-pq), 



qp — pq = 2i sin t R. 



Substituting (|5.68D into (|5.67D and using ( |5.66| ) give 

e**i?(i) = TVpg -sinti?'(i). 

Now 

q =^ — ^ — q + Rq, p = — ^ — p + Rp. 



(5.66) 

(5.67) 
(5.68) 
(5.69) 



Using these and (5.68) we get 
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{2isintRy 



.N+1 



q + Rq] p + q[i 



N -1 



-p + Rp 



.N- 1 



p + Rp ] q — p [ —i 



.N+1 



q + Rq 



iN {qp - pq) . 



This and (5.68) may be written 



and ( 5.69 ) may be written 



{sintky = N^{qp), 
svatR — 3 {qp), 

e'*R{t) + sintR'{t) = Npq. 



Thus (taking the square of the absolute value of both sides) 



{costR{t) + sin tR'{t)f + sm^tR{tf = (sintR) 



I 2 



Using 



jx — —a • JX It — —11, It — — — 



8~R' 



which follow from ( ^!66|) , we get the second order equation 



R{ty + 2 sint costR{t) R! {t) + R' {tf 

1 /I R" (t)^ 
= - I -sin^i +smt costR'Ut) + (cos^t + N'^ sin^t) R' (t) 

2 V4 R'{t) 

From the Neumann expansion of the resolvent kernel R{t] A) — R{t) we obtain the expansion 



(5.70) 



R{t, \)=po + 2plt + Apie + -(1 - iV^ + 36^2)^2^3 ^ 



where 



Po = 



27' 



(5.71) 



If we denote by Eiq{Q\s) the p roba bility that an interval (of the unit circle) of length s contains no eigenvalues 
(modifying here the notation of then 



£;^(0;s) = exp|-2^ R{x,l)dx^ 



(s = 2t). 



Using the expansion (5.71) with A = 1 (and additional terms computed from the differential equation ( 5.70| )) we find 
that 



^Jv(0;^5) 



1 - s 



36 



+ 



17640 



1 
14 



675 



1_ 10 A 



2N^ 2N^ 



s« + • • • (s ^ 0) 



where we have replaced s by so that the N ^ oo limit is clear. This converges uniformly for all N and bounded 
s. Observe that the corrections to the limiting coefficients are 0(1/ N'^) as N ^ oo. 
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VI. GENERALIZATIONS OF HERMITE, LAGUERRE AND JACOBI 



In this final section we shall show that there are differentiation formulas of the form (1.4) for very general orthogonal 
polynomial ensembles, and that if the weight function is the standard Hermite, Laguerre, or Jacobi weight function 
multiplied by the exponential of an arbitrary polynomial then the coefhcients m{x), A(x), B{x) and C{x) in ( |1.4| ) are 
themselves polynomials. Some, but not all, of our derivation can be found in the orthogonal polynomial literature 
but our presentation will be self-contained. 

Throughout, we shall write our weight function as 

w{x) = e-^("). 

As stated in the Introduction, the polynomials orthonormal with respect to w{x) are denoted Pk{x) {k = 0, 1, • • •), 
and we set ipk{x) = pk{x)w{xY^^ so that {ipk} is orthonormal with respect to Lebesgue measure. The underlying 
domain T) of all these functions we take to be a finite or infinite interval. We are interested in differentiation formulas 
(1.4) when, up to constant factors, 

ip{x) ^ iPn{x), i>{x) ^ ipN-i{x). (6.1) 

It is well-known that if fcjv denotes the highest coefficient in pn{x), and if 

kN-i 
«Af = — ; , 

kN 

then there is a recursion formula 

xpn{x) ^ aM+iPN+i{x) + bNPN{x) + aNPN-i{x) (6.2) 
as well as the Christoffel-Darboux formula 

JV-l 



Pk[x)Pk[y) = aN . (6.3) 

x — y 



(See Chap. 10 in ||.) 
We shall always assume that our weight function satisfies 

x^w{x) is bounded for each fc = 0, 1, • • • 

and that V{x) is continuously differentiable in the interior of D. And we assume that 



lim w{x) = (6.4) 

x^dD 



although this wih be relaxed later. We define 



ANix)=aN ipN{y)'^N-iiy)U{x,y)dy, Bn{x) ^ qn I (pNiy)^U{x,y) dy . 
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Proposition: // (6.4) holds and and are given by l[6.y then we have (1-4) with 



m{x) = 1, A{x) = -An{x) - \v'{x), 



B{x) = Bn{x), C{x) 



BN^lix). 



Proof. If 
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iV-1 



then for any polynomial n{x) of degree at most — 1 

7r(x) = / K{x,y)TT{y)w{y)dy. 



Apply this to tt{x) — p'j^{x) and integrate by parts, using (S.4) to eliminate any boundary terms. We obtain 

dK{x,y) 



dy 



PN{y)w{y)dy + / K{x,y)V'{y)pN{x)w{y)dy. 



Now both dK /dy and A', as polynomials in y, are both orthogonal (with respect to w) to PN{y) since they have 
degree at most N —1. It follows that the first integral above vanishes and that we can write the resulting identity as 



p'n[x)^ J K{x,y) {V'[y)-V'{x))pN{y)w[y)dy 

= -aN j [pN{x)pN-i{y) - PN-i{x)pN{y)]pN{y)U{x,y)w{y) dy 
by (p^) and the definition of U{x,y). We have shown 



p'n{^) = ~An{x)pn{x) + Bn{x)pn~i{x). 
It follows from this, of course, that 

PN-ii^) = -An-i{x)pn-i{x) + Bn-i{x)pn~2{x). 
We use the recursion formula (^^) (which holds also if each p^ is replaced hy (Pn) and find that this is equal to 

•fN-iiy) [y(pN-iiy) - anfNiy) - bN-iVN^i{y)\U{x,y) dy pn-i{x) 

+ / fN-iiyW{x,y)dy [xpN -lix) - aNPN{x) - bN~iPN-i{x)] 



(6.5) 



= / y^N-iiy)i^ - y)Uix,y) dy pN-iix) + ANix)pN-iix) - 



BN-lix)pN{x). 



The last integral on the right side equals 



P%-i{y) iV'ix) - V'{y)) w{y) dy = V'ix) + / p%_,{y)w'iy) dy 



The last integral vanishes, as we see by integrating by parts and noting that p'j^^i is orthogonal to pn-i- Thus we 
have shown 



p'n_i{x) = [An{x) + V'{x)]pN^i{x) 



aN 



BN~lix)pNix). 



(6.6) 



The statement of the proposition now fohows from (6.5) and (3.6) if we use the fact that 



fN = Pn'w + t;—Pnw ' , 
2 w 



and similarly for ip'pf_i. 



Remark. The assumption (6.4) is not just a technical requirement. The conclusion of the proposition is false without 
it. (Consider, for example, the Legendre polynomials on (—1, 1), where V{x) — an d the conclusion of the proposition 
reads ip'n{x) = 0.) Nevertheless, we shall be able to handle some cases where (|6.4[) fails. 



Example 1 (generalized Hermite). Here V{x) is a polynomial of even degree (at leas t 2) with positive leading coefficient 
and V — (—00,00). The conclusion of the proposition holds and so we have (L4) with m{x) = 1, with A{x) a 
polynomial of degree at most degV^ — 1, and with B{x) and C{x) polynomials of degree at most deg — 2. 
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Example 2 (generalized Laguerre). Here 



X e 



-W(x) 



where a > — 1 and is a polynomial of degree at least 1 with positive leading coefficient, and T) — (0,oo). In this 
case 



a , W'{x)~W'{y) 

U{x,y) = h 

xy X ^ y 



(6.7) 



Now ( p.4| ) is satisfied if a > and the proposition tells us that in this case we have (1.4) with 

m{x)^x, A{x)^-xAn{x) + ---W'{x), 



B{x)=xBn{x), C{x) 



a-N 



xBn-i{x) 



(6.8) 



now A{x) is a polynomial of degree at most deg W while B{x) and C{x) are polynomials of degree at most deg W —1. 

To extend this to all a > —1 we see that there are problems in the integrals defining Aisf{x) and Bn{x) arising 
from the term a/xy in (6.7). The contribution of this term to the integral defining An{x), say, equals (we assume 
now a > 0) 



ajv- 



PNiy)pN-iiy)e-'^^y'>y"~Uy. 



Integration by parts shows that this equals 



X 



{pN{y)pN-i{yy 



y"dy. 



This expression is well-defined for all a > —1 and in fact represents a function of a which is real-analytic there. (The 

coefficients of the pn are clearly real-analy tic f unctions of a.) 

This argument shows that both sides of (L4), with the coefficient polynomials given by (6^), are (or extend to be) 
real-analytic for a > — 1. Since they agree for a > they must also agree for a > —1. 



Example 3 (generalized Jacobi). Here 

w{x) = (l-x)"(l + x)'^e~'^(^) 
where a,/3 > —1 and is a polynomial, V = (—1, 1). In this case 

a(l +x){l + y)+ (3{1 - x){l - y) W'{x) - W'{y) 



U{x,y) = 



(1 - a;2)(i _ 



and the proposition tells us that for a,/? > we have (1.4) with 

m{x) ^l-x\ A{x) = -(1 - x^)An{x) - |(1 + x) + ^(1 - x) - l^W'{x), 
B{x) = (1 - x')Bn{x), C{x) = ^^(1 - x^)Bn-i{x)- 

now A{x) is a polynomial of degree at most deg -I- 1 while B{x) and C{x) are polynomials of degree at most deg W. 
The identity can be extended to a, > — 1 as in Example 2. (It is convenient to express the integrals defining Aiq{x) 
and Bn{x) as sums of integrals by using a representation 

1 = u{x) + (1 - u{x)) 

where m = 1 in a neighborhood of a; = — 1 and u = in a neighborhood of x = 1; this separates the difficulties at the 
two end-points. The details are left to the reader.) 

Remark. It is clear that the last examples can be generalized to any weight function of the form 
where W{x) is a polynomial and for each which is in the closure of V we have > —1. 
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